Fluctuations of the correlation dimension at metal-insulator transitions 
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We investigate numerically the inverse participation ratio, P2, of the 3D Anderson model and of 
the power-law random banded matrix (PRBM) model at criticality. We found that the variance of 
lnP2 scales with system size L as o^{L) — (t^(oo) — AL~^'^^^'^ , being D2 the correlation dimension 
and d the system dimension. Therefore the concept of a correlation dimension is well defined in the 
two models considered. The 3D Anderson transition and the PRBM transition for h = 0.3 (see the 
text for the definition of h) are fairly similar with respect to all critical magnitudes studied. 



PACS number(s); 71.30.+h, 72.15.Rn, 73.20.Jc 

Critical eigenfunctions at metal-insulator transitions 
show multifractality, whose study has been very inten- 
sive in the last two decades [0,^. Each critical eigen- 
function (/>Q(r) can be characterized by a set of inverse 
participation ratios (IPR) 



j |0„(r)p« d'^r, 



(1) 



where the index a label the different eigenfunctions and 
d is the embedding dimension. In a good metal the IPR 
scale with size L as F,^"^ oc while in an insulator 

Pq"^ OC i". Wegner found, from a renormalization- 
group treatment of the metal- insulator transition in 2 -t- e 
dimensions, that the average IPR at criticality show an 
anomalous scaling of the form 

Pq cx L--D<j(«-i) , (2) 

being Dq a set of generalized fractal dimensions. 

The IPR fluctuations were studied for 2D systems in 
the framework of the supersymmetry method |^-|^ and 
it was found that the distribution function of Pq°''^ nor- 
malized to its typical value p^^'^'^ is scale invariant at 
criticality and inversely proportional to the squared adi- 
mensional conductance [Q. Although the 2D case does 
not present a true Anderson transition, the previous re- 
sult motivated the conjecture that in general the distri- 
bution function of Pq"''^ normalized to its typical value, is 
universal, i.e., size independent for L 00. Accordingly, 
it is assumed that the distribution function of InPg""* 
is a universal curve which is just horizontally shifted by 
changes in L. 

Recently it was claimed by Parshin and Schober ||^ , on 
the basis of numerical calculations, that the correlation 
dimension D2 in the 3D Anderson model at criticality is 
not well-defined due to strong fluctuations in the IPR. 
For the sizes considered, it seemed that the standard de- 
viation of the distribution of In P2"'' grows with system 
size proportionally to InL. Then, for L — > 00 the corre- 
lation dimension, instead of tending to a single value, it 
would tend to a universal distribution. 



If confirmed, the previous result would impose dras- 
tic changes in our understanding of critical properties 
at metal-insulator transitions. Mirlin and Evers 
addressed this problem and study theoretically and nu- 
merically the fractal properties of the power-law random 
banded matrix (PRBM) model at criticality. They found 
the distribution function of the IPR to be scale indepen- 
dent for all the values of the parameter b characterizing 
this model (see below). The model describes a whole 
family of critical theories parameterized by 5, in the same 
way as the dimensionality labels the different Anderson 
transitions. The standard 3D Anderson transition should 
be equivalent to a PRBM model with a 6 of the order 
of unity, although this belief is not based on any direct 
knowledge. They finally claimed that the disagreement 
between their results and those of Rcf. ||] is due to the 
small system sizes used in ^ and not to the different 
models employed. 

Our aim is to perform a careful statistical analysis on 
the data from numerical calculations of both the 3D An- 
derson model and the PRBM model at criticality using 
system sizes larger than in previous calculations [^,0 . We 
want to study the system size dependence of the fluctua- 
tions of the IPR and to elucidate: i) whether the correla- 
tion dimension at the Anderson transition is well-defined 
or alternatively presents an scale invariant distribution, 
ii) whether this transition is equivalent or not to a PRBM 
model at criticality. 

We first consider the standard Anderson model on a 
3D simple cubic lattice, represented by a tight-binding 
Hamiltonian with matrix elements of the form 



(3) 



where i,j denote lattice sites. The diagonal (site) ener- 
gies are randomly distributed with constant probability 
in the interval —W/2 < < W/2, and the off-diagonal 
elements tij are taken equal to unity for nearest neigh- 
bors, which sets the energy scale, and to zero otherwise. 

The PRBM model, introduced in Ref. ||^, describes a 
ID sample with random long-range hopping. It can ap- 
proximately represent a variety of physical systems from 
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an integrable billiard with a Coulomb scattering cen- 
ter |10[ , to the Luttinger liquid at finite temperatures 
p] , |l2| . The model is represented by real symmetric ma- 
trices whose entries are randomly drawn from a normal 
distribution with zero mean and a variance depending on 
the distance of the matrix element from the diagonal 
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This model was shown to undergo a sharp transition at 
a = 1 from localized states for a > 1 to delocalized states 
for a < 1. This transition is supposed to be similar to 
an Anderson metal-insulator transition, presenting multi- 
fractality of eigenfunctions and non-trivial spectral com- 
pressibility at criticality. The parameter b determines the 
critical dimensionless conductance and so establishes the 
character of the transition. For 6=1 the nearest level 
spacing distribution differs from the typical one at the 3D 
Anderson metal-insulator transition Recent calcu- 

lations by us have obtained different diverging exponent 
for the correlation length as the transition is approached 
from below and from above [ p^ . 

We obtain the eigenfunctions and eigenvalues of the 
Hamiltonian matrix by numerical diagonalization. In the 
case of the Anderson model, we use techniques for large 
sparse matrices, in particular a Lanczos tridiagonaliza- 
tion without reorthogonalization method | p5[ , while for 
the PRBM case we employ standard diagonalization sub- 
routines, since we have to deal with full matrices. For the 
Anderson model, the system size varies between 5 and 40, 
and we consider a small energy window (—1, 1) around 
the center of the band. We take for the critical disorder 
the value Wc — 16.5. In the PRBM case, the system size 
ranges between L — 100 and 15000 and the energy win- 
dow considered is (—0.4, 0.4). Reducing the width of the 
previous windows do not alter the results. The number 
of random realizations is such that the number of states 
included for each L is roughly equal to 3 x 10^, except 
for L = 40 in the Anderson model where this number is 
2 X 10**. In order to reduce edge effects, we use periodic 
boundary conditions in all cases considered. 

To elucidate whether the correlation dimension D2 
possesses a well defined single value or alternatively cor- 
responds to a distribution, we have calculated the IPR 
for the wavefunctions in the energy window considered 
for many disorder realizations. For each L we obtain the 
distribution function of lnP2"\ since this is a well be- 
haved self-averaging magnitude. At the end of the paper, 
we will also discuss the distribution of In Fig. 1 we 

show the evolution of the distribution function F(ln Pj"' ) 
with L for the 3D Anderson model. The system sizes 
drawn are L = 5, 6, 8, 10, 12, 15, 17, 20, 25, 30 and 40, 
from right to left. It is clear that ^(InPj"'') changes with 
size. As L increases this distribution becomes wider and 
the height of its peak smaller. This is in qualitative, but 
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FIG. 1. Distribution function F(lnP; 
son model on a logarithmic scale for L 
17, 20, 25, 30 and 40, from right to left. 



for the 3D Ander- 
5, 6, 8, 10, 12, 15, 



not quantitative, agreement with the results of Ref. 
We characterize the previous distribution by its av- 



erage value. 



and its variance, cr^(L) 



var(lnP2 ). This variance increases with L and seems 
to saturate at a constant value, as one can implicitly ap- 
preciate from the peak of the distributions. We try a fit 
of the form 



(5) 



with (T^(oo), A and 7 being three adjustable parameters. 
We found that the exponent 7 was always very close to 
the correlation dimension divided by 2d. Thus, we fix 



= £1 
^ 2d 



(6) 



and keep only two free parameters. In the language of 
scaling theory, this exponent characterizes the behavior 
of the irrelevant length. Our assumption that the irrele- 
vant exponent 7 is equal to D2/2d properly interpolates 
between two known limiting cases of the PRBM model. 
For & <C 1, one can get from Eqs. (||) and (^) the b\nL 
correction predicted in Similarly, for 6 ^ 1, the expo- 
nent 7 in Eq. tends to 1/2, which implicitly coincides 
with the results of Ref. [|| for this regime. 

In Fig. 2 we represent on a log- log scale {00) — a^{L) 
as a function of L for the 3D Anderson transition. The 
fitted values of the free parameters are ct^(oo) = 1.09 
and A — 1.24. The slope of the straight line has not 
been fitted and corresponds to D2/2d, where D2 = 1.4 
has been obtained from the inset of Fig. 2. In this in- 
set we plot (luPj"-*) versus InL. The straight line is a 
linear fit to the data and its slope is equal to D2. The 
value cr^(cxD) = 1.09 found is in good agreement with the 
conjecture a'^{oo) « 1 |3|. 



2 





I 




-1.1 



t> -1 3 



-1.5 



-1.7 




-2.0 



-2.5 



-3.0 



-3.5 



InL 



10 



-4.0 



InL 



FIG. 2. (j^(oo) — o'^{L) as a function of L for the Anderson 
transition on a log-log scale. L ranges between 5 and 40. The 
straight line is a linear fit to Eq. (^) , and the slope is obtained 
from the inset. Inset: {InPg""*) versus InL. The straight line 
is a linear fit to the data and its slope is equal to the (average) 
correlation dimension D-z- 

These results are in contradiction with those of Parshin 
and Schober who found a hnear increase of (j{L) with 
InL. We have checked that (t(L) versus InL do not fol- 
low a linear behavior. We believe that the disagreement 
is due to two reasons. Firstly, the relatively small sys- 
tem sizes employed in Ref. and secondly the use of 
different estimates of an effective IPR, instead of directly 
studying the width of the InR^"^ distribution. On the 
other hand, Mirlin and Evers H already obtained for the 
PRBM model similar results to ours. 

Eq. (||) may be also valid for other transitions, like the 
integer quantum Hall transition, where the irrelevant ex- 
ponent y -which is the same as our 7- is known to be 
y = 0.4 ±0.1 and D2 — 1.48. These values are consistent 
with Eq. d^). Polyakov ||l^ and Evers and Brenig |l^] 
already obtained a relation between the irrelevant expo- 
nent y and the correlation dimension D2, but our result 
seems to fit more naturally the reported values of y. 

Although D2 is well defined in the macroscopic limit, 
it is clear from Fig. 2 that the results of numerical 
calculations with present computers will drastically de- 
pend on the particular definition of correlation dimen- 
sion adopted. For the sizes available the distributions 
of lnP2^"^ are so wide that one can obtain very differ- 
ent numerical values of D2 depending on the method of 
calculation, as pointed out in Ref. but this does not 
imply at all the existence of a distribution of correlation 
dimensions in the macroscopic limit. 

We have now performed similar calculations for the 
PRBM model, considering several values of the parame- 
ter b in the range where the transition could be similar 



FIG. 3. 0-^(00) — o-'^(L) as a function of L on a log- log 
scale for the PRBM transition with 6 = 0.1 (squares) and 0.3 
(diamonds). Solid lines are fits to Eq. (^). 

to the 3D Anderson transition. As regards the IPR fluc- 
tuations, the results for the PRBM model follow the same 
behavior, given by Eq. (^, as the Anderson transition. In 
order to show this, we plot in Fig. 3 the same quantity as 
in Fig. 2, but for the PERM model with h = 0.1 (squares) 
and with b — 0.3 (diamonds). The data in all cases are 
well fitted by straight lines, whose slopes are equal to 
L'2/2rf. The fitted parameters for the cases shown in 
Fig. 3 are (t^(cx)) = 0.55 and A = 0.56 for b = 0.1, and 

(00) = 0.33 and A = 0.35 for 6 = 0.3. The values oi D2 
have been obtained by the same procedure as for the An- 
derson model and are equal to 0.21 and 0.48 for b = 0.1 
and 0.3, respectively. For 6=1, <j'^{L) is practically con- 
stant, explaining the scale- invariance claim made in Refs. 

For 6 > 1, Eq. (^) is still valid, but A changes sign. 
In other words, now the system approaches the asymp- 
totic value ct^(cxd) from above. 

Now we focus on the behavior of the asymptotic value 
ct^(oo) as a function of 6, which is summarized in Fig. 4. 
We were not able to fit all the data (including large val- 
ues of b not shown in the figure) with a simple function. 
Power laws of 1/6 with different exponents are able to 
fit extended parts of the curve. For example, the dashed 
line corresponds to a 1/6 dependence and fits fairly well 
the intermediate regime, 6 1, in which we are most 
interested. For 4 < 6 < 12, the results are in reasonable 
agreement with the 1 /6^ prediction of Evers and Mirlin 
0, based on the fluctuations of the 2D case and tak- 
ing into account that the dimcnsionless conductance is 
proportional to 6. 

The distribution of 2: = Pj" V-F2*'*'''^ where the typical 
IPR p!^^^^^ is chosen as the median of the distribution 
Ip^ , is very different from the distribution of In P2°'^ due 



3 




FIG. 4. Asymptotic value of <y^{L) as a function of h. The 
dashed line is proportional to 1/6. 

to the their strong asymmetries and long tails. We have 
studied the size dependence of the variance ct^ of z. For 
6 ^ 1 we have analytically calculated this dependence 
using the renormalization group method of Levitov . 
One combines two blocks of size L to form a block of 
size 2L. Some states are assumed to remain practically 
unchanged in this process, while others combine to form 
a more extended state. From the evolution equation de- 
scribing the behavior of the distribution function, one 
can obtain D2 in terms of the parameters of the model 
1^. We extended the procedure to obtain the following 
equation for the variance 



:(2i) 



b 0-2-D2 



(7) 



We have checked that for b < 0.3 Eq. (0) fits the variance 
of z fairly well. On the other hand, the numerical cal- 
culations for larger values of b and for the 3D Anderson 
transition do not follow this behavior. Eq. (Q) predicts 
a power law approach of tT^(i) to its asymptotic value 
with an exponent —D2/S, in contrast with the exponent 
— 02/2(1 observed for the variance of lnP2^"^ 

We have compared the Anderson model with the 
PRBM model for several values of b through different 
critical parameters. The case b = 0.3 is very similar to 
the Anderson transition, presenting practically the same 
values for all the critical magnitudes studied. The cor- 
relation dimension divided by the embedding dimension 
is equal to 1.4/3 = 0.47 for the Anderson model and to 
0.49 for the PRBM model. Correspondingly, the scaling 
with system size of the fluctuations of InPj"-* is rather 
similar. The asymptotic value of these fluctuations for 
the PRBM model is the same as for the Anderson model 
divided by d. Also, the normalized variance of the near- 
est level spacing at criticality 1 20 2^ , which characterizes 
the intermediate statistics at the transition |23l, is 0.18 



for the Anderson model and 0.19 for the PRBM model 
with b = 0.3. 

Following the comparison between the two models, we 
have finally analyzed the large values tail of the distribu- 
tion of z. Mirhn and Evers (7|j| predicted for the PRBM 
model a power-law tail F{z) oc where X2 depends 

on the transition considered. For example, X2 is equal to 
4.2 for 6 = 1 and 2.1 for b = 0.3. We have checked that 
this prediction is also obeyed by the Anderson model, for 
which we obtained a value of X2 equal to 1.6, in relative 
agreement with the PRBM model for b = 0.3. 

In summary, we found that the variance of the fluctua- 
tions of In P2°'^ tends to an asymptotic value, implying a 
well-defined correlation dimension. For the intermediate 
regime of the PRBM model, cr^ (00) is proportional to 1/6. 
For the Anderson transition and the PRBM model at 
criticality, the variance tends to its asymptotic value as a 
power law ij-D2/'2-d_ xhus, the irrelevant length scale di- 
verges with an exponent £)2/2d, a result that may be ex- 
tended to other transitions like the integer quantum Hall 
transition. The behavior of the distributions of lnP2"'' 

and of P2^°V^2*''''' is different. We checked that the 3D 
Anderson transition and the PRBM transition for b = 0.3 
are similar for all the critical magnitudes analyzed: the 
correlation dimension, the scaling of the fluctuations of 
lnP2"'', the exponent characterizing the long z tail of the 
distribution F{z) and the normalized nearest level vari- 
ance. 
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